Heisenberg position-momentum uncertainty relation beyond central potentials 
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Recently it was shown in [New J. Phys. 8, 330 (2006)] that the three dimensional Heisenberg 
uncertainty relation arOp > 3h/2 might be significantly sharpened if the relevant quantum state 
describes the particle in a central potential. I extend that result to the case of states which are 
not the eigenstates of the square of the angular momentum operator. I derive a new lower bound 
for arCTp which involves the mean value and the variance of the square of the angular momentum 
operator. 
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I. INTRODUCTION 

The famous Heisenberg uncertainty relation P, is 
one of the most fundamental results of the whole quan- 
tum meehanies. Furthermore, testing various uneertainty 
relations provides a way to characterize many impor- 
tant quantum features. For a generic three dimensional 
system described in terms of a wave function ^ (r) the 
Heisenberg uncertainty relation reads 



3^1 



(1) 



where ct^ = \/{r'^) - (t) • (r) and cTp = \/{p'^) - {p) ■ (p) 
denote the standard deviations of position r and momen- 
tum p variables respectively. 

In a special case where the wave function is an eigen- 
state of the square of the angular momentum operator 
(r) = ti^l {I + (r) the lower bound in ([T]) was sig- 
nificantly sharpened in [3| (this result has been recently 
rediscovered in Q) 



(p2) >h[l + 



(2) 



The fact that (r) = = {p) for all eigenstates of 
was taken into account. We shall also notice that all 
results presented in Q were derived in a general case of 
Z?-dimensional systems. 

The eigenstates of the operator play a special role 
in atomic physics, where a quantum system is usually as- 
sumed to evolve in a central potential depending only on 
the distance from the origin = rr. For several central 
potentials the product Crfp was analytically derived and 
checked numerically in [5|. Various uncertainty relations 
are of special importance in numerical computa- 
tions related to electronic structures, in particular to that 
based on the density functional theory , because they 
provide sensitive tools to verify physical adequateness of 
obtained electron densities [l5|. However, since the re- 
sults of numerical calculations and possible experiments 



cannot guarantee us that -tp (r) has all expected proper- 
ties, we shall ask about the validity of the uncertainty 
relation ^ when ip (r) is not exactly the eigenstate of 
L^. The aim of this paper is to give an answer to that 
question. 

Obviously, the standard deviations and dp are in- 
variant under translations in both position and momen- 
tum (transformations r i-^ r + tq and p i— >■ p + Po), 
however the average value of is not. Therefore, let me 
introduce an invariant version of the angular momentum 
operator (r and p are operators) 



(3) 



This is the angular momentum in the position and mo- 
mentum reference frames centered at (r) and (p) respec- 
tively. Let me then assume that we have at our disposal 
the average value of the square of this angular momentum 
operator described by a dimensionless parameter Ljnv 



^ d\r{r)Ll,^{r). (4) 



In Section |TT1 I present a simple example which shows 
that when one docs not know the state (one knows 
only the parameter Linv) one cannot refine the general 
Heisenberg bound ([T|). This happens because even if 
Lhiv ~ / (^ + 1) for some / £ N, one can still construct 
states laying far away from the eigenstate of 1/?^^ labeled 
by the quantum number I. To overcome this issue I shall 
employ the variance of L^^^ and, in addition to Ljnv, use 
the dimensionless parameter 



(5) 



Of course TZinv = only for the eigenstates of Lf^^. Thus, 
relatively small values of TZinv shall justify the approxi- 
mation that V' is some eigenstate of Lf^„. 

The main result of this paper is the sharpened version 
of the Heisenberg uncertainty relation ([Ij 
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derived in Section IIIII In the discussion part (Section 
IIV[) I show that this new uncertainty relation hnks in 
a continuous manner the previous result ([2|) recovered 
from ([6]) for Ljnv = \Jl {I + 1) and TZinv = 0, with the 
Heisenberg uncertainty relation ([T|) achieved when Ljnv = 
or T^iiiv — > oo. We shall notice that for (r) = = (p) 
we obtain L?^^ = what provides a full correspondence 
with the case ([2]). 



Heisenberg-Iike uncertainty relation for photons. The 
most important advantage of this method is that one 
does not need to rely on commutation relations between 
the conjugate variables. In the first step I define the fol- 
lowing functionals: 



(11) 



II. THE LOWER BOUND (gl) BEYOND 
EIGENSTATES OF L^^ 

In this Section I show that the constraint on the aver- 
age value Lf^^ (from now on we put h = 1) of the L?^^ 
operator (or simply operator in the reference frame 
where (r) = = (p)) does not lead by itself to an in- 
equality stronger than the Heisenberg uncertainty rela- 
tion ([!]). To this end let me use the function 



fi ir,0,^) 
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7ri/4Y {21 + 1)1 



r'e~^'I^Y^{BM. (7) 



which is the normalized eigenstate of some isotropic har- 
monic oscillator with the angular momentum quantum 
number Z, the magnetic number m = 0, and the principal 
quantum number n ~ 0. Now we shall take the superpo- 
sition of the ground state /o (angular momentum number 
equal to 0) with the state for (^o + 1) > Lfnv ^'^'^ 
> 1 



The state 



'^o(Zo + l) 



^0 (^0 + 1) 

satisfies (r) 
but 
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V^o (?o + 1) 
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0, (p) = and (L? 
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This observation means that keeping (^Lf^-^^^ constant we 

arc able to go arbitrarily close to the ground state fo tak- 
ing arbitrarily large Iq- However, for the example state 
([5]) we can check that 

7^inv = Lf„^ (^0 (^0 + 1) - LI„) — > oo, (10) 

which explains, why the Heisenberg lower bound 3/2 can 
be asymptotically reached in ©. 



III. PROOF OF THE UNCERTAINTY 
RELATION ^ BASED ON CALCULUS OF 
VARIATIONS 



P'[r,t^]= d'rr{r){-A)i^{r), (12) 



where A is a three-dimensional Laplacian and we put 



A. Ordinary Heisenberg uncertainty relation 

Following the idea presented in [l^ I shall briefiy de- 
scribe the variational method using the example of the 
Heisenberg uncertainty relation ([IJ. However, I modify 
this derivation and use the Lagrange multiplier to include 
the normalization condition. To prove ([T|) one can start 
with solving the following variational equation 



S 



[t/.*, V'] [V-*, V'] - a / rf'V i^*{r) ^ [r) -1 



0, 



(13) 

where A is the Lagrange multiplier related to the nor- 
malization constraint of the wave function. The equation 
(HH) gives 



A V (r-) = 0, 



(14) 



where we denote u = XP. In the next step we impose 
the normalization condition (we multiply (|14p by ip* and 
integrate over cPr) and find that A = 2ur' . Finally, we in- 
troduce a dimensionlcss variable ^ = r / X and obtain the 
Schrodinger equation for the three dimensional isotropic 
harmonic oscillator 



1 



-A 



(15) 



The eigenvalues of the left hand side of ([15)) are 
a;(n-|-3/2), n G N, thus, the eigenvalue equation for 
(ITSl) gives 



uj = n-\- 



(16) 



The smallest possible value Wmm of uj in ([T6|) is equal to 
3/2 for n = and when ^ (^) is the ground state. This 
result provides us the inequality 



XP > LOr, 



(17) 



A main ingredient of my derivation shall be the vari- 
ational approach used recently in [l6| to prove a new 



Now we choose the coordinate and momentum refer- 
ence frames such that (r) = = (p). In this specific 
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frame wc have X 
reads 



(Tr and P = Op, and our result (H?]) two additional constraints: (L 



> 



(18) 



The observation that (|18|) is invariant under translations 
(is valid in all reference frames) completes the proof 
of ([T]). In other words, in order to prove the Heisen- 
bcrg uncertainty relation ([T]) one needs to show that 
(p2) > 3/2, when (r) = = (p). 



const and 



= TZ ^ const. Note that in this part 

of the derivation I use the ordinary angular momentum 
operator L. The variational equation in that case reads 
(F = 7^ + L^) 



B. Derivation including constraints related to 
angular momentum 

In order to prove the main result of this paper (|6]) I 
shall use the method described in IIII Al together with 
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09) 



and leads to the equation 

(^^r^ - X^A - A + 27]P + 2AL4j V (r) = 0. (20) 

The parameters A, 77 and A in ([T5)) play the role of La- 
grange multipliers and the factor —2 before rj and A was 
introduced for further convenience. Imposing the nor- 
malization constraint we find that A = 2Lu'^ + 2riL'^ + 2AF. 
Thus, we have to solve the following counterpart of the 
Schrodinger equation ([T5|) 

(-^A + ^ + VL' + AL*) ^ (0 = £i' (C) , (21) 

where £ ^ uj'^ + ryL^ + AF. 



1. Solutions to the eigenproblem 

In the first step we notice that all eigenstates $ (^) of 
the isotropic harmonic oscillator [l3| 

(^-i A + $ (I) = ^ + 2n + z) $ (0 . (22) 

are also solutions of (pij) . However, taking V' (4) to be 
one of these eigenstates, we will not be able to fulfill 
the constraints on angular momentum (in particular all 
eigenstates $ (^) have 71 = 0). Thus, in order to solve 
the equation (|2ip together with the constraints we shall 
look for the solution of (PT]) in the form of a superposition 



of the eigenstates of the harmonic oscillator 

N N 

We assume that each two states (f>i and (for i 7^ j) in 
(|23p differ by at least one quantum number, i.e. 7^ rij 
or li 7^ Ij. Thus, all states $i (^) in the superposition 
([^ are independently the eigenstates of (PT|) with the 
energies 

f., = Q + 2n, + +7^k (k + l) + All {k + 1)' . (24) 

In that way we obtain a set of N independent equations 
that must be simultaneously satisfied 

=£, z = l,2,...,iV. (25) 

But, in our problem we have only three constants to spec- 
ify: £, rj, A (it is important to notice that amplitudes 
Ci do not appear in thus, allowed solutions ([25)1 

of ((2T|) can be the superposition of at most three eigen- 
states of the harmonic oscillator. In other words only the 
value = 3 allows us to find the solutions t] {u!,ni,li), 
A{Lo,ni,li) and £ {ui,ni,li). I will not write them down 
explicitly, but I restrict myself to give in IIII CI a simple 
derivation of lo consistent with these solutions. 

It might appear unusual that we have not specified 
the Lagrange multipliers using the related constraints, 
but during the process of solving the variational equa- 
tion ([T^ . However, we still have at our disposal three 
amplitudes Ci, C2 and C3 and we shall use these coeffi- 
cients to fulfill the remaining constraints. 
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2. Solutions to the constraints 



C. Final uncertainty relation 



For the sake of simplicity let me introduce the following 
notation: a = h {h + 1), (3 = I2 {h + 1), 7 = ^3 (^3 + !)■ 
The three constraints we have imposed lead to three 
equations for the probabilities |Ci|^: 



|Ci| 



I Co 



|C,|^ = 1, 



a|Ci|'+/3|C2r+7|C3r = L2, 



a2|Ci|V/?2|C2|V7'|C3|'=F. 
The solutions of (|26all26c() are: 



\C2 



(a - /3) (a - 7) 

2 _ q;7 - (g + 7) + F 
(/3-a)(/3-7) ' 



|^^|2 ^ g/? - (g + /3) + J- 
(7 - a) (7 - /5) 



(26a) 
(26b) 
(26c) 

(27a) 
(27b) 
(27c) 



In that way we have specified the moduli of Ci , C2 and 
C3 coefficients, but their phases might be arbitrary. The 
solutions (|27aj|27cp possess a permutational symmetry (a 
permutation oiU produces the permutation of |Ci|), how- 
ever at this point we shall decide about some hierarchy 
between the U numbers. Let me choose h > h 1^ h what 
means that a > /3 > 7. Since the solutions (|27all27c() 
must be positive we obtain the following conditions for 
a, /3 and 7: 



/37 - (/3 + 7) + > 0, 
a7 - (a + 7) + F < 0, 
a/3 - (a + /3) + F > 0. 



(28a) 
(28b) 
(28c) 



The conditions (|28all28c() might be also presented in 
terms of two, mutually exclusive cases: 



0<7<L2</3< 



F-7L2 
L2-7 



< a, 



(29) 



or 



0<7</3<L2< 



L2 -7 



F- pi? 
<«<^^. (30) 



In fact, from the beginning we could expect that either 
one or two parameters among (a, /3, 7) shall be greater 
than 1?. 



Since all constraints are fulfilled we are able to de- 
rive in a simple way the coefficient uj consistent with the 
equations (j25p . To this end I write down these equations 
explicitly (for i e {1,2,3}) 

Lo[^ + 2n,+ l}j+7Mk + l)+Alf{h + if ^uj^+f]L''+AF, 

(31) 

multiply i-th equation by |Ci| /cu and sum up over i. 
Due to the conditions (|26aj|26cp all the terms with the 
Lagrange multipliers 77 and A present on both sides of 
pip cancel, and we immediately obtain the solution 



(32) 



The final step is to find values of the quantum num- 
bers rii and It that minimize w. Since the probabilities 
\Ci\^ do not depend on rii we shall take the lowest levels 
Jii = 0, as in the case of the ordinary Hcisenberg uncer- 
tainty relation. An answer to the question which values of 
the parameters li minimize the function uj (li) = uj (0, U) 
is more subtle, because we have assumed that e N. 
However, for a moment we shall treat U as continuous 
parameters (we assume that U € R) and calculate the 
following derivatives: 



|Cir>0, (33a) 
|C2|'<0, (33b) 
|C3|'>0, (33c) 



duj {k) 


ih 


-h){h-h){s~ 


hZi) 


dh 




Q12Q13Q23 




duj (li) 


ih 


-h){h~h){s^ 


'h) 


dh 




Q12Q13Q23 




doj ik) 


ih 


-h){h-~h){S' 


M3) 






Q12Q13Q23 





where S ^ 2 + h + h + h S'lid Qjk ^ 1 + Ij + h- 

We have defined uJmin as a minimal value of w (h) with 
the assumption that h € N. Thus, since w (h) increases 
with h and h and decreases with h we can obtain a lower 
bound Q < LJmin taking the smallest possible values of h 
and ^3 and the largest value of h, according to the ranges 
([29l [30| . This means that the optimal values of G M 
arc: 



h=h 



V4F -f L2 1 
2L 2' 



h 



0, 



(34) 



and lead to the result 

3 T 4 



2 2F(L,7^) 



L2 



(35) 

When we substitute in ([35]) the value F {L,n) ^^TZ + L^ 
we obtain the inequality 



3 L3 VUTiUTl^ - L , , 

XP>UJrmn>n^- + — ^ ^^^^5 . (36) 
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IV. DISCUSSION 




5 10 15 20 



FIG. 1: (Color online) The dependence on Linv of the new 
bound ((6]) for some fixed values of TZim'- Ti-im ~ 10'^ - 
red/dashed, TZ\m ~ 10* - green/dashed-dotted, TZim ~ 10^ 
- blue/dotted. The black line represents the reference bound 
All bounds are above the Heisenberg bound 3/2. 



In this Section I would like to present the connection of 
the new uncertainty relation ^ with the previous results 
([T]) and First of all we shall note, looking at the 

expression (|36|. that 

^^(0,'^) = l> and lim ^2(L,7^) = ^. (41) 

2 7?.— >-oo 2 

These results mean that in both cases, when the angu- 
lar momentum is L = or when the variance TZ is very 
large, we have no improvement of the ordinary Heisen- 
berg uncertainty relation ([1]) . This conclusion is in a full 
agreement with logical expectations. Furthermore, it is 
easy to check that 

1}(V^^TT),0) =^+i, (42) 



Similarly to the result pT)) . the uncertainty relation ([36 
in the reference frame defined by (r) = = (p) reads 



l3 Vl^ 

2 



- 4L4 + 47e - L 



(37) 



what coincides with the uncertainty relation The 
value TZ = Q means that the state is a true eigenstate of 
(or Lf^^ in a general reference frame) and the aver- 
age value must be equal to ^ (Z -I- 1) where I is related 
quantum number. Fig. [T]summarizes these observations. 



where now, since L = iinv, we can replace L by iinv and 
the couple of parameters (L,72,) by (Linv, 7?.inv). Finally 
we shall use the fact that i/jnv is invariant under trans- 
lations both in positions and momenta, thus the uncer- 
tainty relation p7p is also invariant and the main result 
^ of this paper is proven. 

In fact, we are able to improve immediately the in- 
equality ([6]) because we can take an exact, but more 
complicated value Wmm instead of its lower bound f2. 
If U € N, then the optimal values of U arc: = 0, and 



V4~F + I? 
2L 



\/4F + L2 
2L 



(38) 

where [•] and [-J denote the integer valued ceiling and 
floor functions respectively [l^. Finally we have 



2U 



-mv -T Linv) 



(39) 



where the function W (x) reads (i^inv = T^u 

LL ((1 + r^i )' + (! + bJ )' + r^i N - 1) -^^i. 



W(x) 



(i + M)(i + Lxj) 



(40) 

One can check that this result is actually a minor im- 
provement of ([6]), moreover the function W{x) is not even 
continuous. 



V. CONCLUSIONS 

I this paper I have discussed the Heisenberg uncer- 
tainty relation for position and momentum with addi- 
tional information about angular momentum. I derived 
new lower bound for the product of standard deviations 
CTrUp which depends on the average value and the vari- 
ance of the Lf^^ operator. This operator is a square of 
the angular momentum operator defined in the reference 
frame where (r) = = (p). I showed that the relation 
(|n]) links the ordinary Heisenberg uncertainty relation ([T]) 
with stronger relation ([2]) valid for the eigenstates oi L?^^. 
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